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ABSTRACT 

We report the discovery by the Palomar Transient Factory (PTF) of the tran- 
sient source PTFllagg, which is distinguished by three primary characteristics: 
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(1) bright {Rpeak = 18.3 mag), rapidly fading (A_R = 4 mag in At = 2d) optical 
transient emission; (2) a faint {R = 26.2±0.2mag), blue {g'—R = 0.17±0.29mag) 
quiescent optical counterpart; and (3) an associated year-long, scintillating radio 
transient. We argue that these observed properties are inconsistent with any 
known class of Galactic transients (flare stars. X-ray binaries, dwarf novae), and 
instead suggest a cosmological origin. The detection of incoherent radio emis- 
sion at such distances implies a large emitting region, from which we infer the 
presence of relativistic ejecta. The observed properties are all consistent with 
the population of long-duration gamma-ray bursts (GRBs), marking the first 
time such an outburst has been discovered in the distant universe independent 
of a high-energy trigger. We searched for possible high-energy counterparts to 
PTFllagg, but found no evidence for associated prompt emission. We therefore 
consider three possible scenarios to account for a GRB-like afterglow without 
a high-energy counterpart: an "untriggered" GRB (lack of satellite coverage), 
an "orphan" afterglow (viewing-angle effects), and a "dirty fireball" (suppressed 
high-energy emission) . The observed optical and radio light curves appear incon- 
sistent with even the most basic predictions for off-axis afterglow models. The 
simplest explanation, then, is that PTFllagg is a normal, on-axis long-duration 
GRB for which the associated high-energy emission was simply missed. However, 
we have calculated the likelihood of such a serendipitous discovery by PTF and 
find that it is quite small (~ 2.6%). While not definitive, we nontheless specu- 
late that PTFllagg may represent a new, more common (> 4 times the on-axis 
GRB rate at 90% confidence) class of relativistic outbursts lacking associated 
high-energy emission. If so, such sources will be uncovered in large numbers by 
future wide-field optical and radio transient surveys. 

Subject headings: stars: flare - stars: gamma-ray burst: general - stars: super- 
novae 



1. Introduction 

From accreting stellar-mass black holes in our Galaxy to distant active galactic nuclei 
(AGNs) and gamma-ray bursts (GRBs), outflow velocities approaching the speed of light are 
common in nature. Indeed, the number of known sources capable of generating relativistic 
ejecta has expanded in rec ent years to include a c ore-collapse supernova without an accom- 



panying GRB ('SN2009bb: ISoderberg et al.ll2010f). as we l l as the presumed tidal disruption 



of a star by a supermassive black hole (JLevan et al.ll201ll : iBloom et al.ll201ll : IZauderer et al. 
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201ll : iBurrows et al.ll2011t ICenko et al.ll2012l ). With revolutionary new time-domain facili- 
ties slated to come online in the coming decade, even more exotic examples will surely be 
uncovered. 

Time- variable high-energy emission (X-rays and 7-rays) tends to be the hallmark of 
such relativistic outflows. Yet there is good reason to expect that some relativistic outbursts 
may lack a detectable high-energy signature. In the case of GRBs, for example, the most 
mundane possibility is a lack of sky coverage: the most sensitive high-energy GRB detectors 
cover only a fraction of the sky at an y given time. But other, more interesting possibilities 



exist, including viewing-angle effects ( iRhoadsl 119971 : iPerna fc Loeblll998: iNakar et al 



and some physical process suppr essing the high-energy emission entirely (JDermer et al. 



Huang et al. 



2002 



1) 



2OO2I 



2000; 



Rhoadsll2003[ ). The search at longer wavelengths for these "orphan" (i.e.. 



off-axis) afterglows or "dirty fireballs" has remained one of the most sought-after goals in 
the GRB field for more than a decade. 

In this work, we report the discovery by the Palomar Transient Factory of PTFllagg, a 
rapidly fading optical transient associated with a year- long, scintillating radio counterpart. 
The detection of a faint, blue, quiescent optical source at the transient location suggests 
a cosmological origin for the transient (i.e., well beyond the Milky Way and any nearby 
galaxies). At such distances, the observed radio emission requires the presence of relativistic 
ejecta. 

Throughout this work, we adopt a stan dard ACDM cosmol ogy with Hq = 71 km s^^ 
Mpc~^, flm = 0.27, and ^a = 1 — ^m = 0.73 (jSpergel et al.l 120071 ). All quoted uncertainties 
are la (68%) confidence intervals unless otherwise no ted, and UT times are used throughout. 
Reported optical magnitudes are in the AB system (lOke fc Gunnlll983l ). We have corrected 
the reported optical and ne ar-infrared (NIR) ph otometry for a foreground Gala ctic extinction 
of E(B — V^) = 0.044 mag (jSchlegel et al.lll998l ). using the extinction law from lCardelli et al. 
fll989h . 



2. Discovery and Basic Analysis 
2.1. Optical/Near-Infrared 

2.1.1. Observations 

Regular monitoring observations of field 100033 (centered at a = 08*^23^32.42*^, 5 = 
+21°33'34'.'5, with a total on-sky area of 7.2 deg^) were obtained with the Palomar 48 inch 



Oschin telescope (P48) equipped with the refurbished CFIIT12k camera (JRahmer et al. 
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P48 R-band 
2011 Jan 30.22 



« 



Keck/LRIS g-band 
2011 Sep 26.60 * 
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^ 2011 Sep 26.60 
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Fig. 1. — Optical imaging of the field of PTFllagg. The P48 discovery (i?-band) image is shown in 
the left panel. Follow-up Keck/LRIS g-band observations, obtained on 2011 September 26, are displayed 
in the center (wider field) and right (zoomed in) panels. The location of PTFllagg, as determined from 
our P48 imaging, is indicated with a solid circle (1" radius; note that this is significantly larger than the 
astrometric uncertainty in our alignment between the Keck/LRIS and P48 images, which is '^ 50 mas in each 
coordinate). A faint, unresolved source consistent with the location of PTFllagg is detected in both our 
5-band and i?-band (not shown) images. All images are oriented with North facing up and East to the left. 



20081) as part of a program to study stellar variability in Pra esepe (the Beehive Cluster 



Agiieros et al.l l201ll ) by the Palomar Transient Factory (PTF; iLaw et al.l l2009l : iRau et al. 
20091 ). Over 500 individual P48 frames, each with an exposure time of 60s, were obtained 



over the period from 2009 November through 2012 March. All P48 images were obtained 
with a Mould R-hand filter , which is similar to the r' filter f rom the Sloan Di gital Sky Survey 
(SDSS; lAiharaet ahlboill ). but offset by ~ 27 A redward Jofek et al.lboi^ . 



In an image beginning at 5:17:11 on 2011 January 30, we detected a bright but short- 
lived optical flare at the (J2000.0) location a = 08^22"^17.195^ 5 = +21°37'38'.'26, with 
a la astrometric uncertainty of 70 mas in each coordinate (Figure [1]). This source was 
subsequently dubb ed PTFllagg by our automated discovery and classification pipeline 
(JBloom et al.ll2012l ). Our P48 photometry of PTFllagg, calculated with respect to nearby 
point sources from SDSS, is presented in Tabled] 

The peak observed magnitude, obtained in our first image of the field on 2011 January 
30, was measured to be i? = 18.26 ± 0.05 mag. In the next ten P48 images of the field, all 
obtained on 2011 January 30, the source is seen to decay by 1.2 mag in the R band. A faint 
detection is also obtained by coadding all P48 images from 2011 February 1 {R = 22.15±0.33 
mag). The resulting P48 -R-band light curve is plotted in Figured All subsequent P48 images 
result in nondetections at this location. 

Examining our pre-outburst (i.e., before 2011 January 30) P48 imaging, we find no 
evidence for emission at this location in any individual frames (extending back in time to 
2009 November). The typical limiting magnitude for an individual P48 image is i? > 20 mag. 
Stacking all frames from 2011 January 29 (i.e., the day preceding discovery), we limit the 
optical emission at the location of PTFllagg to i? > 21.9 mag. Similarly, coadding all 
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Fig. 2. — Optical light curve of PTFllagg, compa red with a repre sentative sample of afterglows of long- 



duration GRBs discovered by the Swift satellite (C enko et alJ 120091 ). The Swift GRBs are color-coded by 



redshift; small black points indicate GRBs with unknown distance. The observed power-law decline from 
PTFllagg (a = 1.66) is consistent with GRB afterglow observations at Ai « Id after the burst. Though 
at the high end of the observed brightness distribution at At « 0.2 d, a sizeable fraction (~ 10%) of Swift 
events have a comparable i?-band magnitude at At a: 1 d. Inverted triangles mark 3a upper limits. 



pre-outburst P48 images results in a nondetection with R > 23.7 mag. 

Deep optical imaging of the location of PTFllagg was obt ained at late tim es (At > 
1 month) with the Low Resolution Imaging Spectrometer (LRIS; lOke et al.lll995l ) mounted 
on the 10m Keck I telescope {g'- and i?-band filters), and the Inamori-Magellan Areal 
Camera and Spectrograph (IMACS; JDressler et al.ll201ll ) mounted on the 6.5m Magellan- 
Baade telescope at Las Campanas Observatory (J-band filter). 

In our deepest epoch of post-outburst optical imaging (2011 September 26 with Keck/LRIS, 
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or At = 240 d), we identify a faint, unresolved (in 0'.'6 seeing) source in g' and R at (J2000.0) 
coordinates a = 08'^22'^17.202% 6 = +21°37'38'.'26 (Figure [I]). Given the uncertainty in 
the astrometric tie between the Keck/LRIS and P48 imaging (50 mas in each coordinate), 
the observed 90 mas radial offset is not statistically significant (null probability of 0.17). 
Coadding Keck/LRIS images of the field of PTFllagg from several individual nights with 
less ideal conditions (2011 March 4, March 12, and April 27), we can recover an object at this 
location with similar brightness in both g' and R. No emission is detected at this location 
in the /-band IMACS images to / > 25.2 mag. 

We obtained near-infrared (NIR) imaging of the location of PTFllagg with the 1.3 m 
Peters Automated InfraRed Imaging TELescope (PAIRITEL; iBloom et al.l |2006[ ) on 2011 
March 1 (At = 30d). A total exposure time of 2246 s was obtained simultaneously in the 
J, H, and Kg filters. Raw data files were processed using standard NIR reduction methods 
via PAIRITEL Pip eline III (C Klein et ai, in preparation), and resampled using SWarp 
( JBertin et al.l 120021 ) to create I'.'O pixel"^ images for final photometry. 



We also observed the location of PTFllagg with the Wide-Field Infrared Camera 
(WIRC; I Wilson et al.l |2003| ) mounted on the 5 m Hale telescope at Palomar Observatory. 
Images were obtained in the Ks filter on 2012 March 28 (At = 423 d) for a total exposure 
time of 1200 s. The individual frames were reduced using a custom pipeline within the IRAF 
environment ( JTodyl Il986l ) . Both the PAIRITEL and WIRC images were calibrated with 
respe ct to bright field stars from the Two Micron All-Sky Survey (2MASS; ISkrutskie et al. 
2006h . 



No emission was detected at the location of PTFllagg in any of the NIR bandpasses. 
The most constraining limits come from the WIRC observations {Ks > 22.6 mag). 

A full listing of our optical and NIR photometry is presented in Table HJ To convert the 
Ve ga-based measureraents fr om 2MASS to the AB system, we have used the offsets derived 
bv lBlanton fc RoweisI J2OO7I ). 



2.1.2. Constraints on the Decay Index and Optical Outburst Onset 

We fit the observed P48 detections on 2011 January 30 and February 1 to a power-law 
model of the form f^, = /o(i^ — ^o) °; where f^ is the flux density (in /iJy), to is the time of 
the outburst onset, a is the power-law index, and /o is the flux density at a fiducial time 
(to + Is). We find best-fit values of a = 1.66 ± 0.35 and to = 23:34 UT (±1.7hr) on 2011 
January 29. We note that the inferred outburst onset to occurs 16.6 hr after the preceding 
P48 nondetection on 2011 January 29 (Table H]). 



In the event that PTFllagg is a bona fide GRB-hke afterglow (Q, an ahernate con- 
straint on the explosion date can be derived by comparing the peak brightness of PTFllagg 
with the observed distri bution of GRB optical afterglows. Using the comprehensive sample 
from iKann et al.l ( 120101 ). an observed magnitude oi R = 18.26 at discovery implies an age 
of At < 0.5 d (Figure |2]). Put differently, the brightest known GRB optical afterglows reach 
an observed magnitude of i? ~ 18 approximately 12 hours after the onset of the high-energy 
emission. Together with the P48 nondetection on 2011 January 29.31, we can conservatively 
constrain the outburst onset to fall within the window from ~ 17:00 on 2011 January 29 to 
5:17 on 2011 January 30 (55590.71-55591.22 MJD). 

While the overall power-law fit quahty is acceptable (x^ = 8.1 for 9 degrees of freedom), 
we caution that the early opti c al light curves of GRBs rarely exhibit single power-law decays 
JPanaitescu fc Vestrandl bood : iRvkoff et allbood : lOates et al.1l2009[ ). In the event that the 
outburst occurred later than our derived to, the true power-law index will be smaller than 
what we have inferred, and more consistent with most previously observed GRB optical af- 
terglows. If the outburst actually occurred earlier, the decay index would steepen somewhat. 
But temporal indices a > 2.5 are ruled out based on the nondetection on 2011 January 29. 



2.1.3. Likelihood of Quiescent Source Association 



Here we wish to estimate Pchance, the a posteriori likelihood that the coincident quiescent 
counterpart detected at late times in our Keck/LRIS imaging is unrelated to PTFllagg 
(i.e., the transient source). We have measured the areal surface density of objects of this 
brightness in our imaging of field 100033, finding a{R < 26.2) = 0.0 3 galaxies arcsec"^; we 



note that this is consistent with the results from iHogg et al.l 119971 using entirely different 
fields. Using 150 mas, or three times the uncertainty in th e astrometric t ie bet ween the P48 
and Keck/LRIS images, as our search radius, and following iBloom et al.l ( 120021 ). we find that 
-Pchance = 2 X 10~^. Wc therefore consider it highly likely that this source is the quiescent 
counterpart of PTFllagg; however, we consider alternative possibilities below as well. 
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2.2. Radio 

2.2.1. Observations 



We began radio observations of the field of PTFllagg witli the N ational Radio Astr on 
omy Observatory's (NRAO^) Karl G. Jansky Very Large Array (VLA: |Perley et al.ll201ll ) 



on 

2011 March 11 (At = 40 d). The array was in the "B" configuration until 2011 May 6, then 
the "BnA" configuration until 2011 June 1, and the "A" configuration thereafter. Over the 
course of our monitoring, the angular resolution ranged from 0'/3 to 1'.'2. The VLA data were 
reduced with the Astronomical Image Processing System (AIPSJj. For flux calibration, we 
used the source 3C 147, while phase calibration was performed using the objects J0823+2223 
and J0832+1832. As a check of our flux calibration, we have verifled that flux measurements 
of our phase calibration sources remain stable throughout the course of our observations. 

We observed PTFllagg at high frequencies (mm wavelengths) with the Combined Array 
for Research in Millimeter- wave Astronomy (CARMA) beginning on 2011 March 14 (At = 
43 d) and continuing for approximately one month. For our CARMA observations, the 
array was in the "D" conflguration, and the beam had an angular diameter of 10". The 
total bandwidth (lower sideband and upper sideband) was 8 GHz, and the local oscillator 
frequency was 93.6 GHz. The optical depth at high (230 GHz) frequency ranged from fair 
(r ^ 0.4; phase noise ~ 50°) on March 14 and April 11, to good (r ^ 0.1; phase noise ~ 40°) 
on April 7. Data were reduced using standard techniques within the MIRIAD environment 



dSault et al.lll995[ ). 



A transient radio counterpart was detected with both facilities. The radio counterpart 
was unresolved (smallest beam size of 190 mas) and consistent with zero circular polarization 
[q < 10 %) at all epochs. The results of our EVLA and CARMA monitoring are displayed 
in Table |21 while the 8 GHz light curve is plotted in Figure Ei 



2.2.2. Spectral Energy Distribution 

To calculate the radio spectral energy distribution (SED), we must interpolate the var- 
ious observing frequencies to a common epoch. To provide the longest lever arm, we per- 
form this analysis at the two epochs of our 93 GHz CARMA detections: 2011 March 14.05 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation (NSF) 
operated under cooperative agreement by Associated Universities, Inc. 



See http://www.aips.nrao.edu/index.shtml . 
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^— Mean GRB Radio Afterglow 
# PTFllagg at = 1.5 




Time Since Explosion (d) 



Fig. 3. — The 8 GHz radio light curve of PTFllagg, at an assumed redshift of 1.5 (in the middle of our 
allowed range: 0.5 < z < 3.0; §4.ip . For comparison, we have plotted the mean lo ng- duration GRB radi o 
light curve (solid gray line), as well as the 25%-75% distribution (gray shaded region uChandra fc Frailll2012l ). 
The variability superposed on the secular decline is likely due to interstellar scattering by electrons in the 
Milky Way, and is not intrinsic to the source. For comparison, at z = 0.5, the 8 GHz spectral luminosity 
would be a factor of 15 smaller, while at z = 3.0 a factor of 6 larger, than the values plotted here. The 
inverted triangle marks a 3cr upper limit. 



(At pa 43d) and 2011 April 7.03 (At ^ 67 d). We have linearly interpolated flux-density 
measurements made immediately before and after these epochs at frequencies of 5 and 8 GHz. 
Due to the relatively sparse coverage at 22 GHz, we have simply adopted the flux density 
at the closest epoch in time (note that for 2011 March 18 we averaged the two 22 GHz 
measurements obtained on this day). The resulting SEDs are plotted in Figure HI 

We fit a power law of the form f^ = f^v^ to the data, where fy is the flux density (in 
/iJy), z/ is the observing frequency (in GHz), (3 is the power-law spectral index, and /o is 
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Fig. 4. — PTFllagg spectral energy distribution at radio frequencies. The observations at lower frequen- 
cies have been interpolated to common epochs {At sa 43 and 67 d) to match the times of our CARMA 
observations. 

the flux density at a fiducial frequency of 1 GHz. For the first epoch (At ^ 43 d), we find 
P = 0.28 ± 0.08. On the second epoch (At ^ 67d), we measure (3 = 0.46 ± 0.07. Given 
the relatively large degree of variability (see below), together with the sparse coverage at 
high frequencies, we adopt /3 = 1/3 as an approximate spectral slope in the radio for the 
remainder of this work. 



2.2.3. Angular Source Size 

The presence of nonthermal radio emission provides two powerful and independent 
means to constrain the angular size of the emitting region. To begin with, the brightness 
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temperature ( Tp) of an incohe r ent radio emitter ca nnot exceed its equipartition value of 



B,eq 



by 



10^^ K (JReadheadI I1994J : iKulkarni et al.lll998l ). The brightness temperature is given 



Tb — —J 2~r^' ^^^ 

where c is the speed of hght, ks is Boltzmann's constant, z/ is the observing frequency, fi, is 
the observed flux density, and is the angular diameter of the emitting region. Adopting 
Tb ^ 10^^ K in Equation [1] thus implies a lower limit on the angular diameter of the source: 



e > 2.1 



/. 



/iJy 



1/2 



u 



GHz 



/ias. 



(2) 



As can be seen from Equation [2|, the strictest lower limits on the size of the emitting 
region are derived from observations at the lowest frequencies (assuming a power-law spectral 
index /3 < 2). Using our 4.5 GHz observation on 2011 March 13, we find O > 7/ias. Most of 
our early observations at 5 and 8 GHz yield comparable (though slightly less strict) limits. 

Separately, we can constrain the an gular size of t he source from the detection of inter- 
stellar scattering and scintillation (ISS; lRickettlll990l ). To quantify the degree of variation 
induced by the scattering electrons, we calculate the modulation index. 



mp[u) 



ifu) 



(3) 



where V{fu) is the variance of the flux density (with respect to an assumed model), (a^) 
is the average of the square of the individual measurement uncertainties, and (f,,) is the 
average of the flux density. 

We calculated the modulation indices at 5 and 8 GHz, neglecting higher frequencies due 
to the relative lack of observations. We fit the light curves at both frequencies to a power-law 
model of the form /j, = foit — to)~°', finding best-fit temporal indices of osghz = — 0.09±0.13 
(i.e., consistent with no temporal evolution) and agcHz = 0.56 ±0.06. This power-law model 
then forms the reference which we use to calculate the variance at each frequency. In this 
manner, we find mp{5 GHz) = 0.42 and mp{8 GHz) = 0.26. 



We use the Galactic electron density distribution model of ICordes fc Lazid (120021 ) to 
derive the relevant ISS parameters, namely z/q, the transition frequency between the strong 
and weak scattering regimes. For the line of sight to PTFllagg (Galactic coordinates / = 
202.08°, b = 29.2°), we find z/q = 11 GHz. For a point source, the maximum degree of 
modulation {nip = 1) will occur at this transition frequency. It is therefore not unreasonable 
to expect our observations at 5, 8, and (possibly) 22 GHz to suffer from some degree of ISS. 
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For z/q = 11 GHz, our observations at 5 and 8 GHz will be in the strong scattering 
regime (z/ < uq). Furthermore, given the relatively broad bandwidth of our observations 
(Ai//z/ ^ 0.1), we consider only refractive scintill ation. For a p oint source, the modulation 
index in the strong, refractive regime is given by (IWalkerlll998l ) 



mp[iy) 



17/30 



(4) 



For the line of sight to PTFllagg, we therefore expect a significant degree of modulation for 
a point source at our observing frequencies: mp(5GHz) = 0.63, mp(8GHz) = 0.82. 

For an extended source, the observed modulati on will be redu ced by a factor of (9^/0)^^^, 
where B^ is the size of the Fresnel scattering disk (jWalkerlll998l ). 



e. 



VDuo 



1^0 



V 



11/5 



/las. 



(5) 



where D is the effective distance to the scattering screen (D = 0.78 kpc for the line of sight 
to PTFllagg). If we solve for the angular diameter corresponding to the observed degree 
of modulation at each frequency, we find 9(8 GHz) = 10/ias and 0(5 GHz) = 34/ias. We 
therefore conclude that the angular size of the emitting region at At ^ 100 d is ~ 20 /las. 



2.3. High-Energy 



2.3.1. J -ray Limits 



At the time of discovery, three primary high-energy facilities were monitoring the sky 
to search for the p rompt emission from GRBs. The Third Interplanetary Network (IPN; 
Hurley et al.l |2010[ ) is a group of nine satellites sensitive to high-energy emission. When 
multiple satellites detect a GRB, the sky localization can be reconstructed from light travel 
time constraints. The IPN provides essentially continuous all-sky coverage (i.e., 100% duty 
cycle), with a sensitivity to fiuences (lOkeV - 5MeV) of S"^ > 6 x 10~'^ergcm~^ (at 50% 
efficiency; i.e., half of the GRBs with this fluence are too faint to trigger the IPN det ectors) . 
In addition to the IPN, the Gamma- Ray Burst Monito r fGBM: iMeegari et al.ll2009l ) on the 
Fermi sa tellite, and the Bur st Alert Telescope (BAT; iBarthelmy et al.l 120051 ) on the Swift 
satellite (jGehrels et al.l 120041 ) . also regularly discover a large number of GRBs. The GBM 
detects bursts down to a 8keV - IMeV fiuence of 5*^ > 4 x 10~^ergcm~^, but has a field 
of view of 8.8 sr (the area of the sky unocculted by the Earth in the Fermi orbit) and a 
duty cycle of > 80%. Likewise, the Swift BAT has detected events with 15-150 keV fiuences 
as low as 6 X lO^^ergcm"^, but only observes a field of view of 2sr with a duty cycle of 
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~ 90%. We caution that, for all three facilities, the high-energy fluence required to trigger 
the onboard GRB algorithms depends on the duration of the event; therefore, the above 
sensitivity limits should be treated only as approximate. 

We have searched all three facilities for GRB triggers from the direction of PTFllagg 
over the time period from 17:00 2011 January 29 to 5:17 2011 January 30 (i.e., the outburst 
onset window derived in §2.1.2p . No triggers were reported by any facility in the direction of 
PTFllagg during this ~ 12 hr window. We further conducted a search for untriggered events 
in the GEM data in the energy range 10-300 keV on several different time scales (0.256 s, 
0.512s, 1.024s, 2.048s, 4.096s, and 8.192 s)o. No potential high-energy counterparts to 
PTFllagg were found. 

Given the field of view and duty cycle of the GBM and BAT, there is a significant likeli- 
hood that events below the IPN sensitivity threshold would be missed by both instruments. 
For example, for a GRB with fluence above the GBM sensitivity level (but below the IPN 
threshold), the probability of a nondetection from both instruments is as high as ~ 40% 
(assuming a uniform and independent distribution of sky pointings for the two instruments). 
Given the relatively large window of time required to search (i.e., multiple Swift and Fermi 
orbits), we consider a fluence of 5*^ < lO^^ergcm"^ (i.e., twice the all-sky IPN sensitivity) a 
reasonable limit on any high-energy prompt emission associated with PTFllagg. Given the 
extremely weak correlati on between prompt 7-ray fluence and optical afterglow brightness 
( JNysewander et al.ll2009[ ). this limit is consistent with the known properties of GRBs and 
their afterglows. 



2.3.2. X-ray Limits 



To search for an X-ray counte rpart, we obtained ob servations of the location of PTFllagg 
with the X-ray Telescope (XRT; iBurrows et al.ll2005l ) on board the Swift satellite on 2011 
March 13 (At = 4 2 d). Data were reduced using the automated pipeline described by 



Butler fc Kocevskil (120071 ). No X-ray source is detected at the location of PTFllagg at this 
time. Assuming a power-law spectrum with a photon index of F = 2, we derive a 3a upper 
limit on the 0.3-10 keV flux of /x < 8 x lO'^^ergcm^^ g-^. 

Finally, we note that no historical X-ray em ission has been re ported at this location, 
neither in the i? 05*74 T All- Sky Survey (0.1-2.4keV: IVoges et al.lll999l ) nor in any compilations 



■^Two individual GBM detectors with a significance of 4. Oct and 3.8a above background were required for 
a trigger to register in this search. 
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of known Galactic X-ray sources (accessed via the HEASARCcI and SIMBADO databases). 



3. Comparison with Known Galactic Transients 

The combination of (1) a rapidly fading optical transient {AR > 4 mag in At = 2 d) and 
(2) a faint, blue {g' — R = 0.17 ± 0.29 mag), quiescent optical counterpart makes PTFllagg 
unique amongst the thousands of discoveries by PTF to date. Together with (3) the long- 
lived (At ^ 300 d) radio emission, here we attempt to simultaneously account for these three 
distinguishing characteristics. 

In order to understand the nature of the emission from PTFllagg, we must constrain 
its distance. In this section, we first consider a Galactic origin by comparing PTFllagg with 
known classes of Galactic transients. 

Assuming the faint optical source is associated with PTFllagg (i.e., the quiescent coun- 
terpart), the measured color, g' — R = 0.17 ± 0.29 mag, implies a spectral type of ~ F2 
(Teff ~ 7000 K) for a main-sequence star. More conservatively, adopting our 3cr limit on the 
color {g' — R < 1.04 mag), we can rule out single main-sequence stars with Tes ^ 4500 K 
(i.e., cooler than spectral type K4). Given the observed brightness, a main-sequence star 
hotter than K4 would lie at a distance d > 90kpc. This firmly rules out an association with 
the Praesepe cluster {d ~ 175 pc); in fact, only 4 globular clusters are known to exist at 
such large distances in the extreme outer halo of the Milky Way (e.g., AM 1 at d ~ 120 kpc; 
Madore fc Arplll979l ). In addition to the extremely small source densities this far in the halo. 



the inferred lower limit on the radio luminosity at such a distance {I'L^ > 10 ergs" ) is 
three orders of magnitude larger than the most luminous kn own stellar r adio sources (e.g., 
RS CVn binaries, FK Com class stars, and Algol-class stars; ICiidell 120021 ). 



While a posteriori unlikely, it is nonetheless important to consider that the quiescent 
optical source may be unrelated to PTFllagg. Absent color information, an optical nonde- 
tection, even at the depth of our late-time imaging, is not sufficient to rule out a Galactic 
origin. With their smaller effective temperatures, low-mass stars and brown dwarfs (in partic- 
ular ultracool stars, with spectral type later than M7) emit little flux in the optical bandpass. 
Furthermore, ultracool stars are known to exhibit high-amplitude, short timescale (minutes 
to hours) opti cal and radio outbursts that have in the past been mistaken for extragalac 
tic transients ( JBecker et al.l |2004J : iKulkarni fc Raul l2006l : iMahabal et al.l l2012l : iBerger et al. 



''See http://heasarc.gsfc.nasa.gov . 
^See http://simbad.u-strasbg.fr/simbad . 
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2012ah . 

We can use our NIR limits on the quiescent emission at the location of PTFllagg 
to calculate the minimum distance to an ultracool star as a function of spectral type; in 
other words, for each spectral type, any object closer than this "detectability" distance 
would be identified in our NIR imaging. For spectral types later than ~ M4, the strongest 
constraint is provided by our deepest epoch of Ks-hand imaging: Ks > 22.6 mag. Using the 
observed K^ - hand mag nitudes and di s tance (parallax) measurements for ultracool stars from 
Dahn et al.l ( 12002| ) and iPatten et al.l (120061 ) , we fit a low-order polynomial to calculate the 
absolute i^s-band magnitude as a function of spectral type, Mk^{ST), where ST = 5 for M5, 
ST = 12 for L2, etc. We find the scatter about our derived absolute i^g-band magnitude fit 
is ~ 0.30 mag (i.e., 30%) over the range M5-T8. We then convert the observed peak radio 
flux density (/jy,peak ~ 300 /xJy) to a lower limit on the radio luminosity (I'L^) using these 
distance constraints. 

The resulting luminosity limits, as a function of spectral type, are plotted in Figure |5l 
For comparison, we have also plotted all radio observations of ultracool stars from the lit- 
erature (see figure caption for references). Our luminosity limits are typically at least two 
orders of magnitude larger than the most luminous known ultracool stellar flares. Even 
comparing with the recently detected flare from the T6.5 dwarf 2MASS J1047 +21, by far 
the coolest brown dwarf detected at radio frequencies (JRoute fc Wolszczanl 120121 ). our limits 



require a radio luminosity a factor of > 20 times larger. We furthermore see n o evidence for 



a high degree of circ ular polarization (common to many, though not all, flares; iBergerl 12006 



Hallinan et al.l 120071 ). and the radio emission from PTFllagg is much more long-lived than 
these low-mass stellar outbursts (durations typically of only hours). 

In addition to stellar flares, binary systems where one member is a compact object 
(white dwarf, neutron star, or black hole) are known sources of optical and radio outbursts 
in the Milky Way. Such a system could circumvent two issues with Galactic transients 
we previously identified. First, the energy release d uring the accretion pro cess is more than 
sufficient to power the observed radio flux; Cyg X-3 ( IGeldzahler et al.lll983[ ). for example, has 



reached peak radio luminosities in excess of 10 erg s~ . Second, the presence of an accretion 
disc can alter the optical color of such systems. Accordingly, our previous inference that the 
quiescent counterpart must lie at c? > 90 kpc would no longer be valid. 

We consider first X-ray binaries, where the degenerate pri mary is a neutron star or blac k 
hole. Of particular interest are the subclass of microquasars ( iMirabel fc Rodrigueg|l999l ). 
whose powerful radio jets exhibit apparent superluminal motion (and thus imply a relativistic 
outflow). Due to the lack of a bright quiescent optical counterpart, we consider only low- 
mass systems where the accretion occurs via Roche-lobe overflow from the nondegenerate 
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Fig. 5. — Lower limits on the radio luminosity of PTFllagg (solid black line). For each spectral type, we 
calculate a minimum "detectability" distance using limits from our NIR imaging (i.e., any source more nearby 
would have been detected). We then convert this distance to a lower limit on the radio luminosity based 
on the observed pe ak flux from PTFl l agK. Shown for comparison are radio observations of ultracool stars 



Berger et al 



2001 



200S, 



2012: McLean et al 



from the literature (IBerger et alj 2010 : 
l2008bial . l20ol 



Bergeij|2006t lAntonova et al. 



Burgasser fc Putmanll2005 : 



2007 



Hallinan et al 



2007 : 



Audard et al.ll2007l: lBergeiJl2002l: 



Route fc Wolszczan 



20121 ). The inferred luminosity is several orders of magnitude larger than that of any 



previously observed low-mass star or brown dwarf, either in a quiescent or flaring state. 



secondary (low-mass X-ray binaries, or LMXBs). 

Black hole LMXBs are typically characterized by we ll-defined "states" : correlation s 



between X-ray spectra, X-ray fiux, and radio emission (JRemillard fc McClintockl 120061 ) . 
Radio emission is observed in a well-d efined region of this hardness-intensity phase space 
(jFender et al.l I2004J : iFalcke et al.l l2004l ). In the so-called low- hard state, thou ght to cor- 
respond to low (< O.OlLEdd), radiatively inefficient accretion ( lEsin et al.l 119971 ). relatively 
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steady radio emission from a jet is observed in most black hole X-ray binary systems. Like 
PTFllagg, the radio spectrum is fiat or inverted {fi, oc z/^, with (3 ~ 0-0.5), and circu- 
larly unpolarized. However, a reasona bly tight correlation exists between t he radio and 



X-ray luminosity in the low-hard state (I Corbel et al.l 120031 : iGallo et al.ll2003[). of the form 
/;y(radio) oc /jy(X — ray)°-^. Using the derived formulation from iGallo et al.l fcOOSl ) and the 
observed radio flux, we would expect an X-ray flux of /x ~ 2 x 10~^^ ergcm"^ s~^ (note that 
this estimate is entirely independent of the distance to the source). This is more than an 
order of magnitude above our derived X-ray limits. We further note that while neutron star 
X-ray binaries do not obey the same radio-X-ray correla tion in the hard s tate, the ratio of 
X-ray to radio luminosity is even larger in these sources (JMuno et al.ll2005l ). 



Alternatively, the most luminous radio flares from LMXBs arise as the system transitions 
through the intermediate state into a bright, quasi- thermal outburst ( jet emission at the 
highest X-ray fluxes appears to be largely suppressed; iFender et al.ll2004l ). Unlike the steady 
radio jets in the low-hard phase, this state transition in the accretion f low (from radiat ively 
inefficient, advection-dominated to geometrically thin, optically thick: lEsin et al.l 119971) can 



sometimes cause the ejection of relativistic material (JMirabel fc Rodriguez 



LMX Bs in this state do not always follow the same radio-X-ray correlation (JGallo et al. 



mm. While 



20031 ). the radio spectrum from this extended emission becomes optically thin. The X-ray 
and optical fluxes can rise by several orders of magnitude on a time scale of only a few days 
during these "X-ray novae , " but typically both take many months to return to quiescence 
JTanaka fc ShibazakilE996l : Icharles fc Co3l2006h . 



PTFllagg differs from these X-ray novae in several major respects. Most importantly, 
to reach the intermediate (and, ultimately, high) state where the radio flare is launched, the 
compact pr imary must be a. ccreting material at a substantial fraction of the Eddington limit 
(> 1-10%; lEsin et al.lll997l ). As a result, these outbursts have been discovered exclusively 
by wide-field X-ray or 7-ray satellites. But for any reasonable Galactic distance scale {d < 
lOkpc), our X-ray limits rule out emission at the level of 10~^ LEdd (for a 1 M0 black hole 
or neutron star). While our X-ray observations were obtained 42 d after the initial optical 
outburst, this is comparable to the e-folding time of these systems. As it requires ~ 1 month 
for the disc naass t o accrete onto the neutron star or black hole (the viscous time scale; 
King fc Ritterlll998l ). this time delay alone cannot account for the many orders of magnitude 
gap between our limits and the required X-ray luminosity. 

In addition to the lack of bright X-ray emission, we note several more characteristics 
that distinguish PTFllagg from known X-ray nova outbursts: (1) the radio emission at 
late times remains unresolved, which is difficult to reconcile with relativistic ejecta in our 
Galaxy; (2) the inverted radio spectrum is inconsistent with the optically thin emission 
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expected at this time; (3) the time scale of the optical decay (At < 2 d) is significantly 
shorter than what is observed in X-ray novae (r ^ 20-40 d); and (4) the location, well 
off the Galactic plane (/ = 202°, b = +29°), is inconsistent with the known population of 
LMXBs ( Ivan Paradijs &: Whitelll995l : IWhite fc van Paradijslll996l ). whi ch have a scale height 



of (iz 5-, 1 kpc (although several prominent counterexamples are known; 



Zurita et al.ll2000l : iHvnes et al.ll2000l : lUemura et al.ll2000l : Mirabel et al. 
200ll 



Tomsick et al. 



2001 



1999 



Wagner et al. 



Finally, we consider systems with a white dwarf accreting material from a (low-mass) 
stellar companion or another white dwarf. Short-timescale optical outbursts have been ob- 
served from AM CVn systems (a white dwarf accre ting H-poor material in a short-period or - 
bit), but never with amplitudes larger than 4 mag (JLevitan et al.ll201ll : iRamsay et al.ll2012l ). 
Dwarf novae, t he analogous phenomenon to X-ray novae in LMXBs (e.g., disc instabilities; 



Cannizzdll993l ). have also in the past been mistaken for extragalactic transients ( jRau et al. 



20071 ). In fact, the amplitude and duration of the optical outburst from PTFllagg are not 
unlik e the most extrem e dwarf novae. But dwarf novae rarely exhibit coincident radio emis- 
sion (IBenz et al.lll996l). and the few known ex amples have outburst durations of only ~ 1-2 
weeks ( jBenz et al.lll983t iKording et al.l 120081 ). where the radio emission closely follows the 
optical evolution. 

To summarize, we have examined a variety of known classes of Galactic transients 
(stellar flares, LMXBs, dwarf novae), and found that PTFllagg does not flt neatly into any 
of these categories. It should go without saying that it is entirely possible that PTFllagg 
represents a new type of Galactic outburst, characterized by (1) bright, rapidly fading optical 
emission; (2) a long-lived radio transient; and (3) an extremely subluminous (Mr ^ 11 mag 
for d = 10 kpc) quiescent optical counterpart. Given the broad agreement between our 
observations and the properties of long-duration GRB afterglows ( §5.ip . we do not further 
explore this possibility here. 



4. An Extragalactic Origin: Implications and Comparisons 

Having rejected a Galactic origin for PTFllagg, we now consider the possibility that 
it resides instead at a cosmological distance (i.e., well beyond the Local Group and into the 
Hubble flow). Assuming the quiescent counterpart is indeed the host galaxy of PTFllagg, 
we constrain possible redshifts in §4.11 Even with these crude constraints, the angular size 
derived in §2.2.31 requires the presence of a relativistic outflow ( §4.21) . Finally, we briefly 
compare the observed properties with those of known extragalactic sources capable of gen- 
erating relativistic ejecta in §4.31 and quickly settle on a long-duration GRB-like outburst 
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(i.e., the core collapse of a massive star) as the most plausible explanation. 



4.1. Redshift Constraints 

Assuming the quiescent optical source is related to PTFllagg ( §2.1.3p . we can place an 
upper limit on its distance based on the absence of redshifted H I absorption along the line of 
sight (i.e., the Lyman break). Our g' detection implies that redshifted Lya (Arest = 1216 A) 
falls at an observed wavelength of X-Lya ^ 4800 A (i.e., the middle of the g' filter bandpass). 
This results in an upper limit on the host-galaxy redshift of 2; < 3. 



Alternatively, assumi ng a modest rest-fr ame UV luminosity for the host galaxy (Mi 



uv 



< 



— 16 mag, or L > O.OIL*; iReddy et al.l 120081 ) . we place a lower limit on the host redshift of 



2; > 0.5. A similar lower limit is derived if we compa re the observed -R-ban d brightness with 



that of known host galaxies of long-duration GRBs ( jJakobsson et al.ll2012l ). 



We therefore conclude that the redshift of PTFllagg should fall somewhere in the range 
0.5 < z < 3.0. 



4.2. Evidence for Relativistic Ejecta 

In §2.2.3[ we derived two independent constraints on the angular diameter of the emit- 
ting region from our radio observations: O > 7 /las at Atobs ~ 42 d, and O ^ 20 /ias at 
Atobs ~ 100 d. To convert these to constraints on the outflow velocity, we use the redshift 
limits derived above: 0.5 ^ -z < 3.0 (corresponding to angular-diameter distances of 1.3- 
1.8 Gpc for a concordance ACDM cosmology). Assuming ballistic (i.e., constant velocity) 
expansion, the angular diameter is then given by 

rpct 

dA{l + z) 

where T is the outflow Lorentz factor (F = (1 — /3^)^^/^), c is the speed of light, t is the 
time since outburst (in the observer frame), cIa is the angular-diameter distance, and z is 
the source redshift. 

At z = 0.5, where our limits on the outflow velocity are the weakest, we find F > 1.2 at 
Atobs ~ 42 d, and F ;^ 1.3 at Atobs ~ 100 d. These limits vary httle over our redshift range 
of interest, due primarily to the limited evolution of the angular-diameter distance over this 
range: at z = 3.0, we find F > 1.6 at Atobs ~ 42 d, and F ^ 1.6 at Atobs ~ 100 d. 
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We therefore conclude that, even at this late time, the ejecta powering the transient 
emission from PTFllag g are at least transrelativi stic. For any more realistic form for the 



ejecta deceleration (e.g. jBlandford &: McKeelll976l ). we infer that PTFllagg was initially at 



least a modestly relativistic explosion. 



4.3. Comparison with Known Relativistic Sources 

Only a handful of extragalactic sources are known to produ ce relativistic ejecta: G RBs, 
with initial Lorentz factors as least as large a s several hundred ( iLithwick &: Sarill200ll ). and 
possibly greater than 1000 (lAbdo et al.ll2009[): AGN s. in particular the subclass of blazars, 
with Lorentz factors as large as 50 (iLister et al.ll2009l): and the recently discovered relativistic 



tidal disruption flar e s fTDFs: iLevan et al. 



Burrows et al. 



2012 



201 



Berger et al. 



Cenko et al. 



2011 



Bloom et al.l 1201 It IZauderer et al.l 



20121). with in itial Lorentz factors ~ 10 (JMetzger et al. 



2011 



2012bl : Ivan Velzen et al.l 1201 ll ) . Though with only one or two examples 



to date, the known relativistic TDFs do not appear to vary in the optical on time scales as 
short as those of PTFllagg, where 5t <^ Id. Furthermore, the SEDs of these sources are 
dominated by the soft X-ray (~ 1-10 keV) bandpass, with peak isotropic luminosities as large 
as Lx ~ lO^'^ergs"^. Even at 2; = 3, our X-ray limits ( §2.3.21) imply Lx < 6 x lO^^ergs"^. 

Blazars, however, are known to vary in the optical on short time sc ales {5t < Id), and 
have in the past been mist aken for optically discovered GRB afterglows (IVanden Berk et al. 
2OO2I : iGal-Yam et al.ll2002l ). But signiflcant {mp > 10 %) interstellar sci ntillation is observed 
only very rarely in blazars (< 1% of the population; iLovell et al.ll2008l ). More importantly, 
the degree of optical variability observed from PTFllagg, in particular the amplitude from 
peak to quiescence {/S. R > 8 mag in At ^ 1 m onth), makes this source unlikely to belong to 
any known AGN class ( iMacLeod et al.ll2012l ). 



On the other hand, a long-duration GRB can naturally accommodate all of the observed 
properties of PTFllagg. We flnd that the standard GRB afterglow flreball model can ac- 
curately reproduce the observed optical and radio light curves ( §5.ip . The small initial size 
of the ejecta explains the observed interstellar scintillation, though this should be quenched 
as the blast wave expands relativistically (usually on a time scale of weeks to months). And 
the faint, blue quiescent optical counterpart is consistent with the long-duration GRB host- 
galaxy brightness distribution for z > 0.5 ( jjakobsson et al.ll2012l ). We therefore conclude 
that the most likely explanation for PTFllagg is a long-duration GRB-like (i.e., massive 
star core collapse) explosion, the flrst time such an event has been discovered at cosmologi- 
cal distances absent a high-energy trigger. 
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5. PTFllagg as a GRB: Untriggered, Orphan, or Dirty Fireball? 

Broadly speaking, there are three reasons why a distant, relativistic outburst may lack 
detected prompt high-energy emission. The null hypothesis is a lack of sky coverage (i.e., 
an "untriggered" GRB), as the more sensitive high-energy satellites {Swift and Fermi) have 
only a ~ 60% combined likelihood of detecting any given event ( §2.3.ip . The limiting 7- 
ray fiuence from the only all-sky satellite available (the IPN) corresponds to an isotropic 
7-ray energy release of E^^iso = (2-200) xlO^° erg from z = 0.5 - 3.0. These values are 



not sufficiently low to rule out typical cosmological long-duration GRBs (JButler et al.ll2007l ). 
let alone the class of sublu minous (e.g., GRB 980425 / SN 1998bw-like) events uncovered in 
relatively nearby galaxies (ISoderberg et al.l l2006bl : ICobb et al.l l2006l : iGuetta fc Delia Valle 
20071 ). Without any additional information, the simplest explanation is that PTFllagg is 
an otherwise normal but untriggered long-duration GRB. 

There exist other, more intriguing, possibilities, however. The second possible explana- 
tion for a GRB-like explosion abse nt any high-energy signature is a viewing-angle effect. Due 
to their high degree of collimation ( ISari et al.lll999l : lRhoadslll999l ). the prompt emission from 
most GRBs is beamed away from our line of sight. However, the long-lived afterglow emission 
may nonetheless be visible, either if the region generating the afterglow is less bea med than 
the 7-ray emitting material (i.e., an on-axis orphan afterglow; iNakar fc Pirarul2003l ). or if, as 
expected, the outflow spreads laterally at la te times and i lluminates an increasing fraction 
of th e sky (i.e., an off-axis orphan afterglow; lRhoadslll997l : iPerna &: Loeblll998l : iNakar et al. 



2OO2I ) . The discovery of a bona fide orphan afterglow would provide robust constraints on the 
GRB beaming fraction, still a large source of uncertainty in calculations of the true energy 
release and the all-sky rate of GRBs. 

Finally, a source may lack detectable high-energy emission altogether, either because 
no high-energy photons were produced, or such emission may be unable to escape to d istant 
obser vers due to some internal suppression mechanism. It has long been noted (e.g., iPiran 
2004 . and references therein) that the baryon composition of the relativistic jet in the fireball 
model must be very finely tuned in order to generate any detectable prompt high-energy 
emission (the so-called "baryon loading problem"). Without any baryons in the ejecta, 
the internal shocks thought to power the prompt emission will not forno- But with too 
large a baryon fraction, the jet will not accelerate to a sufficiently high in itial Lorentz factor 



(Fn ^, 20), inhibiting a ny high-energy emission via e -e"^ pair production (JHuang et al.ll2002 



Ghirlanda et al.ll2012l ). Such explosions, dubbed "dirty fireballs," have long been predicted 



^An alternative possibility is t hat the prompt emission is ge nerated by magnetic dissipation in a Poynting- 
flux-dominated outflow (see, e.g.. lLvutikov fc Blandfordll2003l ). 
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( JDermer et al.ll2000t iHuang et al.ll2002l : lRhoadsll2003[ ) to occur as a result of a modest baryon 



loading of the jet; a proton content as small as M > 10"'* Mq will lower the initial Lorentz 
factor sufficiently (Fq ~ Ey:p^/Mc^), yet can still produce the observed broad-band afterglow. 
Distinguishing between a source that produces no high-energy emission whatsoever and one 
in which these photons are unable to escape is clearly challenging - for the remainder of this 
work we shall refer to such objects generically as dirty fireballs or afterglows lacking prompt 
high-energy emission. 

Here we attempt to discriminate between these competing hypotheses through two dif- 
ferent means. First, we distinguish between on-axis and off-axis models by comparing the 
observed optical and radio emission with analytic and numerical predictions for GRB after- 
glow light curves in the fireball model ( §5.ip . Second, we calculate the rate of PTFllagg-like 
outbursts to determine if it is consistent with the all-sky (on-axis) GRB event rate ( §5.2p . 



5.1. PTFllagg and the Fireball Model 



In the standard GRB afterglow fireball model (see, e.g., iPiranI |2004J for a review), rel- 
ativistic ejecta with (kinetic) energy ii^KE sweep up material in the circumburst medium, 
forming a coUisionless shock and accelerating electrons to a power-law distribution of en- 
ergies with exponent p and minimum Lorentz factor '-y^. It is assumed that a constant 
fraction of the total post-shock energy density is partitioned to the electrons (eg) and the 
magnetic field (cb)- These accelerated electrons then emit synchrotron radiation, powering 
the long-lived X-ray, optical, and radio afterglow. 

The observed afterglow spectrum depends on the relative ordering of three critical fre- 
quencies: the frequency where self-absorption becomes important (z^a), the characteristic 
frequency of the emission (z/m), and the frequency above which electrons are able to cool effi- 
ciently through radiation (z/c). We shall assume that all our observations occur in the "slow" 
cooling regime {v^ < Uc), and that the self-absorption frequency falls below the frequency 
range probed by our observations (z^a < 10^ Hz). 

The light curve produced by such emission depends on the radial profile of the circum- 
burst medium into which the shock is expanding. The simplest circumburst medium to 
consider is one in which the d ensity is constan t (p oc r°). This scenario is also referred to as 
an interstellar medium (ISM; ISari et al.l 119981 ). and is parametrized in terms of the particle 



number density no, where p = mpUo gem ^. 

Long- duration GRBs, however, have been conclusively linked to the deaths of massive 
stars (e.g., IWoosley fc BloomI |2006| ). In the late stages of evolution, massive Wolf-Rayet 
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stars are stripped of their outer H and (possibly) He envelopes in a wind, leaving behind 
a signature p oc r~^ den sity profile that shou ld be discernible in the afterglow light curve. 



Wind-like environments (jChevalier &: Lill2000l ) are parametrized in terms of A*, where p 



5 X 10^^y4^,r ^gcm ^. 

Finally, we note that the hydrodynamical evolution also depend s on the geometry of the 



outflo w. GRBs are now widely believed to be aspherical explosions (lRhoadslll999l : ISari et al. 



19991 ). biconical jets with half-opening angle Oy At early times, the jet emission is coUimated 
into a narrow cone (^eff ~ T"^ ^ ^j, where F is the Lorentz factor of the expanding shock) 
due to relativistic beaming. As the shock slows, however, simple analytic solutions suggest 
that lateral spreading of the jet becomes important, and on-axis observers eventually "miss" 
emission from wider angles. This hydrodynamic transition manifests itself as an achromatic 
steepening in the afterglow light curve (the "jet break"), with an expected post-break decay 
proportional to t~^. While more recent numerical simulations have s u ggested a more complex 



picture of the jet-bre ak phenomenon (jZhang fc MacFadyenI 12009 



van Eerten et al. 2010 



Granot fc Piranll2012l ). the assumption of a general light-curve steepening around the time 
when F = 1/6'j remains largely valid. 

At At ^ 0.2 d (the approximate time of discovery), we expect the optical bandpass 
to fall below the cooling frequency (i.e., z/opt < Vc)- If we also assume that these early 
optical data occur before any jet break (see below), the observed temporal decay index 
(ttopt = 1-66 ± 0.35) can be translated directly into the electron spectral index p. For a 
constant-density medium, we find p = 3.21 ± 0.47, while for a wind-like environment, we 
infer p = 2.55 ± 0.46. Electron spectral indices derived from previous obs ervations of long- 
duration GRBs (jShen et al.ll2006l : IStarling et al.ll2008l : ICurran et al.ll2010l ) fall in the range 
~ 2-3, so the large uncertainty makes it difficult to distinguish between the competing 
density profiles solely on this basis. 

While the radio emission is relatively variable at At > 40 d, the approximate radio 
spectral index, /3radio ~ 0.3, implies that the peak synchrotron frequency u^a is not well 
below the radio at this time (or else we would expect /3radio ~ — !)• Conservatively, we 
assume I'mi^t = 40d) > 10 GHz, and fu^{At = 40d) > 300yuJy. For a wind-like medium, 
z/jn oc t~^/^ and f„^ oc t^^/^. Extrapolating back to the time of optical discovery, we conclude 
Um{At = 0.2d) > 3 X 10^3 Hz, and U^{At = 0.2d) > 4 x 10^ pjj. For u^ < u < u^, 
fu oc z/(i^p)/2 ^ z/^*^'^^. Thus, we find that the inferred optical (i?-band) fiux at discovery, 
fv ~ 500 /iJy, is a factor of ~ 3 larger than our observations at this time. For a constant- 
density environment, the peak fiux is constant in time, and so a similar analysis yields a 
self-consistent result. We therefore do not consider a wind-like medium any further. 



Using these general constraints, we have used the software described by Ivan Eerten et al. 
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( 120121 ) to fit the observed optical and radio light curves to afterglow models calculated from 
high-resolution two-dimensional relativistic hydrodynamical jet simulations. In all cases, we 
have assumed a const ant- density circumstellar medium and adopted a fiducial redshift of 1. 
We find that a relatively wide set of parameters is able to reproduce the observationqj, largely 



consis te nt with values der ived from previous GRB afterglow modeling (jPanaitescu fc Kumar 



2001a| ]bl: lYost et al.ll2003l ). although with a somewhat smaller circumburst density (no ^0.1 



cm~^). The best-fit model, assuming the observer is oriented directly along the jet axis (i.e., 
6'obs = 0), is plotted in Figure E] (-Eke = 3 x 10^^ erg, 6-^ = 0.50rad, Uq = 1 x 10"^cm~^, 
eg = 0.1, en =.0-li and p = 2.9). Repeating a similar analysis with the software described 



by lYost et al.l (120031 ) yields qualitatively similar results. The general agreement between the 
fireball model and our optical and radio observations supports our conclusion that PTFllagg 
is most likely a distant, relativistic explosion like other long-duration GRBs. 

One concern regarding the z = 1 models is the implied angular size of the source (B). 
In the best-fit on-axis model, we find that the outflow has an angular size ^ 40 mas at 
At = 40 d. This is somewhat larger than the value inferred from our scintillation analysis, 
but within a factor of two. However, many of the other z = 1 models that provide a 
reasonable fit to the data are likely to be too spatially extended to scintillate strongly at 
At ^ 100 d. One potential solution may be a more distant origin; due to cosmological time 
dilation, an observer-frame time of At = 100 d would correspond to only 25 rest-frame days 
post-explosion at 2; = 3, half the expansion time as inferred at z = 1. Given the large spread 
in acceptable models, however, we do not explore this possibility further here. 

Next we consider limits on the opening angle and observer orientation from the observed 
optical and radio emission. After the jet break, the peak synchrotron flux declines linearly 
with time (in a constant-density environment). Thus, if the jet break occurred well before 
the first radio observations, the large radio flux would be difficult to reconcile with our early 
optical observations. We consider it likely, then, that tj > 40 d. These conclusions are largely 
confirmed by our numerical models, where we find that the opening angle is only weakly 
constrained to be 6j > 0.15 rad. 

In addition to cases where the observer is oriented directly along the jet axis (i.e., 
6'obs = 0), we also have considered more general geometries, where the observer may be 
oriented off-axis, either within (i.e., 6'obs < Oj) or outside (6'obs > 6'j) the jet opening angle. 
For simplicity, we consider only "top-hat" jet geometries, where the jet Lorentz factor is 
given by a step function. For ^obs > ^j, observers will see rising emission until approximately 



^In all cases the predicted X-ray flux is well below the XRT limit (i i2.3.2|) . so we have not included this 
point in our fitting. 
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Time Since Explosion (d) 



Fig. 6. — The solid line shows the best-fit afterglow model (|van Eerten et al.ll2012f) when the rclativistic 
jet is oriented directly along the line of sight to the observer (i.e., 0obs = 0). The dashed curve displays 
the best-fit model when the viewing angle is allowed to vary freely. Given the relatively sparse dataset (in 
particular the lack of X-ray observations), a wide variety of models are able to reproduce the observed optical 
and radio emission. However, we find it impossible to reproduce the observed emission when the viewing 
angle is outside the cone of the jet (i.e., 6'obs > ^j)- We therefore consider it unlikely that viewing angle alone 
can account for the lack of high-energy emission from PTFllagg. 



the time of the jet break, after which the decay will resemble the on-axis case. For observers 
off-axis but within the jet opening; angle, the modifications to the on-axis afterglow light 



curve s will be more subtle ( JGranot et al.ll2002l : IZhang fc MacFadyerul2009l : Ivan Eerten et al. 
2010h . 



The arguments used above to infer tj > 40 d necessarily require that the observer cannot 
be well outside the jet opening angle (or else we would expect to see post jet-break decay). 
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This result is borne out by our numerical modeling, where geometries with 6'obs > 6'j are 
unable to accurately reproduce the observed light curves. Allowing the observer orientation 
to vary as a free parameter, the best-fit afterglow model is plotted as a dashed line in Figure E] 
{Eke = 9 X lO^^grg, 0. = 0.20 rad, ^obs = 0.19 rad, uq = 1 x 10-3cm-^ e^ = 0.04, e^ = 0.2, 
and p = 3.0). 

Finally, we can derive a lower limit on the distance to PTFllagg based only on the 
observed radio evolution. The radio spectrum at At = 67 d, f^ oc z/^, with (3 ^ 0.3, is in- 
consistent with Sedov- Taylor blast-wave evolution. In other words, the outgoing shock wave 
has not transitioned to nonrela tivistic expansion at this point in time. The nonrelativistic 
transition will occur at a time (IWygoda et al.ll201ll ) 



Even neglecting our previous finding of a low circumburst density, we infer a sizeable lower 
limit on the blast- wave kinetic energy: Ek_e ^ 10^*^ erg. Integrating over the observed 8 GHz 
radio light curve, we measure a fluence of S'rad = 2 x 10~^° erg cm~^. For a typical cosmological 
distance {z = 1, or d^ = 2 x 10^^ cm), this corresponds to a radiated energy of -Erad = 
3x 10^^ erg, a typical radiative efficiency for a GRB. But for a Galactic outburst, the radiated 
energy would be many orders of magnitude smaller (ii^rad = 6 x 10'^^ erg at d = 10 kpc). Unless 
the radiative efficiency was incredibly small, we once again conclude that PTFllagg must 
lie at a cosmological distance. 



5.2. The Rate of PTFllagg-like Events 

Our objective in this section is to estimate the number of GRB optical afterglows dis- 
covered by chance (i.e., not as a result of deliberate follow-up observations of a high-energy 
trigger) by PTF. If the likelihood of chance detection of an untriggered afterglow with PTF 
is significant, we will conclude that the rate of PTFllagg-like events is consistent with the 
rate of normal (i.e., on-axis) long-duration GRBs. If this probability is small, then we can 
use these calculations to place lower limits on the observed frequency of PTFllagg-like out- 
bursts (in units of the GRB rate). Given the relatively complex nature of PTF scheduling 



(JLaw et al.l 120091 ). we have conducted a series of Monte Carlo simulations to this end. 



PTF began full operations on about 2009 April 1. We have retrieved a listing of all 
images obtained beginning at this time through 2012 December 31, or over a period of 45 
months. We removed fields at Galactic latitude |6| < 20° (due to the large foreground 
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extinctionjj. The resulting sample includes 129,206 pointings, each covering an area of 
7.2 deg^. The sample comprises 1940 unique fields, each imaged an average of 67 times. 

Since launch, the Swift BATlj detects GRBs at a rate of ~ 90 yr^^. The field-of-view 
of the BAT is ~ 2 sr, and the instrument has a duty cycle of ~ 90%. Thus, the all-sky rate 
for events at the BAT threshold is ~ 630 yr~^. Over the 3.75 yr period of interest, the total 
number of all-sky GRBs is ~ 2360. We note that this is an upper limit to the long-duration 
GRB rate, as we have included short-duration GRBs in this sample as well. 

For each trial, we create a mock catalog of 2360 GRBs. Each GRB is randomly assigned 
a trigger time to (uniformly distributed between 2009 April 1 and 2012 December 31) and 
spatial coordinates a, 6 (isotropically distributed on the sky). To estimate the duration over 
which the optical afterglow would be detectable by PTF, we utilize the sampl e of 29 long- 



durat ion afterglows from the Palomar 60 inch (P60) Swift afterglow catalog (jCenko et al. 



20091 ). These events were selected solely on the basis of visibility to Palomar Observatory, 
so they should represent an unbiased sample of the Swift afterglow brightness distribution. 
For each event in the P60-Swift sample, we have calculated the amount of time following 
the high-energy trigger that the afterglow is brighter than R = 20 mag. These values range 
from < 204s (GRB 050721) to 1.2 d (GRB050820A). Each mock GRB is randomly assigned 
one of the 29 actual "visibility windows" from this sampla^°l. 

For each mock GRB, we then determine if the event occurred within the 7.2 deg^ foot- 
print of any individual PTF image, and, if so, if the time of observation occurred within 
the necessary window during which the afterglow was brighter than 20 mag. The number of 
afterglows detected in each trial (A^grb), together with the number of individual frames on 
which each detected afterglow was brighter than the P48 sensitivity limit (A^Det), were then 
recorded. The results of 1000 individual runs (i.e., different randomly selected groups of 
2360 GRBs) constitute a sufficiently large sample to evaluate the likelihood of serendipitous 
detection of long-duration GRB afterglows with PTF. 

In the 1000 trials conducted, at least one GRB afterglow was detected (i.e., A'grb > 1) 
in 970 instances. Thus, the probability of detecting at least one on-axis afterglow over the 
course of the first two years of PTF is quite high, P(A'grb > 1) = 97%. The expectation 



^Givcn that the primary objective of PTF is the discovery of extragalactic transients, this represents less 
than 10% of the total number of observations. 

^See lhttp://swift. gsfc.nasa.gov/docs/swift/archive/grb_tablel . 

^°For GRBs without any detected optical afterglow (e.g., "dark" bursts), we use the earliest non-detection 
below our sensitivity threshold for the visibility window. If anything, this would bias us to over-estimate the 
expected number of untriggered GRB afterglow detections by PTF. 
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value for the number of afterglows detected is A = 3.3. The distribution of the number 
of afterglows detected in our 1000 trials is reasonably well described by Poisson statistics 
(Figure [7]). In this respect, then, PTFllagg appears to be consistent with a normal on-axis 
GRB. 

However, the field in which PTFllagg was identified (the Beehive cluster) is atypical 
amongst PTF pointings. Most fields are only observed 2 or 3 times per night (multiple images 
are used to identify Solar-System objects). But the Beehive cluster is a "high-cadence" field, 
observed many times (> 10) per night during its observing season. Instead of calculating 
the rate of afterglow detections over the entire survey (i.e., A^grb), a more appropriate 
comparison would limit the scope to similar high-cadence fields. 

We therefore consider on how many individual images each of the 3340 "detected" GRB 
afterglows (in our 1000 trials) were above the P48 limiting magnitude (i.e., iV^et)- This is 
illustrated in Table [31 The vast majority of the afterglows are detected on only one or two 
images (87%). In fact, in our 1000 trials, an optical afterglow was detected on at least ten 
individual images only 11 times (i.e., P(A^Det > 10) = 2.6%). PTFllagg was detected 11 
times on 2011 January 30 with R < 20 mag. 

We can understand this result analytically in the following manner. In the case where 
the integration time {5t) is much smaller than the period over which a transient is visible 
(r), the number of detectable events at any given time will be 

9 = -. 8 

47r, 

where Vt is the field of view (in steradian) and A^ is the all-sky event rate. For long-duration 
GRB optical afterglows, serendipitous detection by PTF will be dominated by the ~ 10% 
of events that remain brighter than i? < 20 mag for r f^ 1 d (certainly this is true for those 
afterglows with A/'oet > 3). Thus, for the PTF project, VL/At: = 1.7 x lO'Sr (7.2 deg^), and 
adopting A/" ~ 0.1 x 630yr~^ and r ^ 2.7 x 10"^ (Id), we find g ^ 3.0 x 10~^ events per 
field. 

The expected number of detected events. A, will then be gA^obs, where A^obs is the number 
of (independent) measurement epochs. Over the two-year period of interest, the number of 
individual P48 images obtained is A^obs(all) = 1.3 x 10^. Thus, we predict A ~ 3.7, in good 
agreement with the results of our Monte Carlo simulations. 

Conversely, we can calculate the relative frequency of high-cadence (A^obs[> 10]) ob- 
servations in our two-year PTF sample by measuring how often each field was observed on 
a nightly basis. The results of this analysis are shown in the far-right column of Table |3l 
As is evident, high-cadence observations with A'"obs(> 10) (i.e., more than 10 observations 
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Fig. 7. — Normalized histogram of the number of serendipitous detections of normal on-axis GRB afterglows 
by PTF in our 1000 Monte Carlo trials. The distribution is reasonably well described by a Poisson function 
with A — 3.3 (solid black line). 



of a field obtained in a single night) occur with a frequency of 1% when compared with 
regular-cadence fields (A^obs[l] + A^obs[2]). 

From this analysis, we conclude that the rate of PTFllagg-like events is inconsistent 
with the rate of long-duration GRBs with 97.4% confidence. Admittedly a number of as- 
sumptions went into this analysis, and one should always be careful with results drawn from 
such an a posteriori analysis. But independent of the exact likelihood, we conclude that the 
probability of untriggered afterglow detection in a high-cadence PTF field is small. Either 
we have been quite lucky, or we may have uncovered a new, more common class of distant, 
relativistic outbursts lacking entirely in high-energy emission. 
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it is crucial to verify, however, that our inferred rate does not violate any other limits on 
short-timescale transients, either from PTF itself, or from previous optical and radio surveys. 
As highlighted above, low-cadence fields are observed significantly more frequently with PTF 
than high-cadence fields like the Beehive. Thus, any short-timescale (At < 1 d) transient 
should be detected in many more iVobs(l) and iVobs(2) fields than high-cadence fields. In the 
case of PTFllagg, repeating the above Monte Carlo simulations for a transient population 
with five times the GRB event rate (but the same optical brightness distribution), we find 
an expected number of detected sources of A = 16.7 in all fields. At first glance, the fact that 
we have not discovered such a population of sources would seem to favor the untriggered 
GRB scenario. 

Here it is important to distinguish between transient detection, by which we mean a 
source is above the P48 sensitivity limit on a given image, and discover y, where a transien t 
is flagged as astrophysically interesting (by software or human beings; iBloom et al.ll2012l ). 
Because of the large number of uncatalogued asteroids near the PTF limit, our software 
requires at least two detections at a given location to flag a source as a bona fide transient 
(e.g., to "discover" the source). Thus, any PTFllagg-like outburst with only a single de- 
tection (A^Det = 1) will never be discovered by our survey. Likewise, there may be subtle 
biases limiting our capability to identify and/or conduct follow-up observations of similar 
short-timescale transients with only a few detections. 

Whether these discovery biases are sufficient to account for the lack of similar sources in 
our low-cadence fields with PTF remains to be seen. We have attempted to search through all 
PTF discoveries that were detected only on a single night (independent of A'Det)LJ, but have 
yet to uncover any additional viable candidates. Ultimately, future wide-field, high-cadence 
optical surveys may be required to resolve this issue. 



Finally, we co mpare our derived rate of PTFllagg-like events with previous searches 



for or 



2006 



Dhan optical (IVanden Berk et a. 



20081 ) and radio ( iLevinson et al 



2002; 



2002 



Becker et a 



Berger et al. 



2004 



2003 



Rvkoff et al. 



2005 



Gal- Yam et al.ll2006 



Rau et al. 



Soderberg et al 



2006aJ) afterglows, to verify that our results are consistent with these limits. T he tightest con- 



strain ts on the rate of relativistic outbursts come from radio surveys, where iGal-Yam et al. 



( 120061 ) derive a limit on the all-sky volumetric rate of GRB-li ke explosions of A^ < 10^ eve nts 
Gpc~^ yr ^- Even assuming an all-sky GRB rate as large ( Guetta fc Delia Valid 2007 ) as 



100Gpc-3yr-i 



Butler et al. 



more recent estimates suggest a significantly smaller value 
20101 ). a population of PTFllagg-like events occurring at a rate of ~ 5 times that of normal 



^^We cannot avoid the requirement of at least two detections, however. Otherwise we would be completely 
swamped with asteroids, which are detected at a rate of thousands per night. 
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GRBs is consistent with these results. Our derived rate is therefore orders of magnitude lo wer 
than the all-sky rate of Type Ibc supernovae {N = 2.6 x 10^ Gpc~^yr~^: iLi et al.l 20 111). It 



may approach the rate of low-luminosity GRBs (ISoderberg et al.l l2006bl : ICobb et al. 



2006 



Guetta fc Delia Vallell2007l ). although this depends both on the assumed beaming correction 



and the true GRB rate. 



6. Summary and Conclusions 

To summarize our results, we report here the discovery of PTFllagg, a rapidly fading 
optical transient with a long-lived, scintillating radio counterpart. Together with the ob- 
served optical and radio light curves, the detection of a faint, blue quiescent counterpart at 
the location of PTFllagg indicates that the transient likely originated in the distant uni- 
verse. Using our measurements of the source size derived from the radio observations, we 
infer that PTFllagg must be powered by a relativistic outflow. These properties are all 
consistent with the population of long-duration GRB afterglows, marking the first time such 
an event has been discovered at cosmological distances without a high-energy trigger. 

Searching various high-energy satellites, we find no potential 7-ray counterpart for 
PTFllagg. We therefore consider three possible explanations that can simultaneously ac- 
count for a GRB-like explosion without any associated prompt high-energy emission: an 
untriggered GRB, an orphan afterglow, and a dirty fireball. 

Using the all-sky rate of GRBs discovered by the Swift satellite, together with a mea- 
surement of their observed optical brightness distribution, we have calculated the likelihood 
of serendipitous untriggered GRB afterglow detection by PTF (April 2009 - December 2012). 
Surprisingly, we found that the a posteriori probability of untriggered GRB afterglow de- 
tection in a high-cadence field like the one where PTFllagg was found (11 observations on 
a single night) is only 2.6%. While we cannot rule out entirely our null hypothesis that 
PTFllagg is an untriggered GRB, this probability is sufficiently low that we consider alter- 
native interpretations as well. 

The afterglow emission from an orphan GRB will rise in flux at early times, as more 
and more of the jet becomes visible due to relativistic beaming effects. Using both analytic 
and numerical formulations, we are unable to reproduce the observed PTFllagg light curves 
unless the observer viewing angle is within the opening angle of the jet. While these models 
assume a relatively simple jet structure, the requirement of rising afterglow emission at early 
times is a robust prediction for all off-axis models. 

A more intriguing possibility is that PTFllagg may represent a new class of relativistic 
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outbursts w ith little or no corresp onding high-energy emission. In much the same way that 
SN2009bb (ISoderberg et al.ll201Cll ) demonstrated that the more nearby, subluminous class of 



GRBs may generate relativistic ejecta yet lack high-energy emission, PTFllagg may play an 
analogous role for the more energetic, cosmologically distant sample of long-duration GRBs. 
Dirty fireballs (i.e., a baryon- loaded jet) are one possible explanation, though alternative 
possibilities surely exist as well. 

In this picture, the inferred rate of PTFllagg-like events must be four times higher 
(90% confidence) than the rate of on-axis long-duration GRBs. When combined with tradi- 
tional core-collapse supernovae and long-duration GRBs, these objects would enable a more 
complete census of the deaths of massive stars, and also provide a probe of the location 
of massive-star formation in distant galaxies without the need for a high-energy satellite 
trigger. 

Regardless of its ultimate origin, we expect such sources to be discovered in large num- 
bers by ongoing and future wide-field, high-cad ence optical surveys s uch as the Catalin a 
Real-Time Transient Survey forake et allboogl). PTF . Pan-STARRS JXaiser et al.lboioh . 
and the Large Synoptic Survey Telescope (jivezic et al.ll2008l ). Furthermore, the discovery of 
PTFllagg bodes well for optical surveys in the future era of gravitational wave astronomy, as 
the electromagnetic counterparts of gravitational wave sources should exhibit largely similar 
observati onal signatures (though they are also expecte d to be associated with more nearby 
galaxies; iNakar &: Piranll201ll : iMetzger fc Bergerll2012l ). 
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Table 1. Optical/Near-Infrared Observations of PTFllagg 



Date 


Telescope/Instrument 


Filter 


Exposure Time 


Magnitude 


(MJD) 






(s) 




55590.30519 


P48/CFHT12k 


R 


540 


> 21.9 


55591.22026 


P48/CFHT12k 


R 


60 


18.26 ±0.05 


55591.22245 


P48/CFHT12k 


R 


60 


18.25 ±0.04 


55591.23391 


P48/CFHT12k 


R 


60 


18.36 ±0.05 


55591.25326 


P48/CFHT12k 


R 


60 


18.51 ±0.08 


55591.26691 


P48/CFHT12k 


R 


60 


18.51 ±0.04 


55591.26800 


P48/CFHT12k 


R 


60 


18.61 ±0.06 


55591.33081 


P48/CFHT12k 


R 


60 


18.53 ±0.17 


55591.36188 


P48/CFHT12k 


R 


60 


18.96 ±0.28 


55591.40604 


P48/CFHT12k 


R 


60 


19.36 ±0.10 


55591.42439 


P48/CFHT12k 


R 


60 


19.46 ± 0.09 


55591.43978 


P48/CFHT12k 


R 


60 


19.51 ±0.10 


55593.40775 


P48/CFHT12k 


R 


420 


22.15 ±0.33 


55594.23819 


P48/CFHT12k 


R 


300 


>21.2 


55621.19100 


PAIRITEL 


H 


2246 


>20.4 


55621.19100 


PAIRITEL 


J 


2246 


>20.6 


55621.19100 


PAIRITEL 


Ks 


2246 


> 19.7 


55624.49 - 55678.28 


Keck I/LRIS 


9' 


6680 


26.63 ±0.33 


55624.49 - 55678.28 


Keck I/LRIS 


R 


5700 


26.28 ±0.28 


55830.60259 


Keck I/LRIS 


3' 


2100 


26.34 ±0.19 


55830.59849 


Keck I/LRIS 


R 


2160 


26.17 ±0.22 


55944.22461 


Magellan/IMACS 


I 


2400 


> 25.2 


56014.27324 


P200/WIRC 


Ks 


1200 


> 22.6 
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Table 2. Radio Observations of PTFllagg 



Date 


Observatory 


Frequency 


Integration Time 


Flux Density 


(2011 UT) 




(GHz) 


(min) 


(MJy) 


Mar. 11.27 


VLA 


8.46 


15.5 


300 ± 23 


Mar. 13.12 


VLA 


4.50 


15.5 


217 ±34 


Mar. 13.12 


VLA 


7.92 


15.5 


375 ± 28 


Mar. 14.05 


CARMA 


93.5 


492.0 


450 ± 140 


Mar. 15.08 


VLA 


4.50 


15.5 


183 ± 36 


Mar. 15.08 


VLA 


7.92 


15.5 


224 ± 30 


Mar. 18.08 


VLA 


4.50 


15.5 


58 ±37 


Mar. 18.08 


VLA 


7.92 


15.5 


171 ±31 


Mar. 18.10 


VLA 


22.46 


10.8 


237 ± 79 


Mar. 18.15 


VLA 


4.50 


15.5 


93 ±35 


Mar. 18.15 


VLA 


7.92 


15.5 


215 ±29 


Mar. 18.17 


VLA 


22.46 


10.8 


460 ± 77 


Mar. 26.22 


VLA 


8.46 


15.5 


277 ± 24 


Apr. 5.17 


VLA 


21.8 


30.1 


271 ± 18 


Apr. 7.03 


CARMA 


93.5 


378.0 


480 ± 120 


Apr. 10.11 


VLA 


4.8 


13.7 


127 ±21 


Apr. 10.11 


VLA 


7.4 


13.7 


155 ± 16 


Apr. 11.01 


CARMA 


93.5 


342.0 


-40 ± 150 


Apr. 18.15 


VLA 


4.8 


37.2 


81 ±11 


Apr. 18.15 


VLA 


7.4 


37.2 


165 ±8 


May 6.08 


VLA 


4.8 


13.7 


232 ± 18 


May 6.08 


VLA 


7.4 


13.7 


221 ± 14 


May 13.15 


VLA 


4.8 


12.9 


117 ±18 


May 13.15 


VLA 


7.4 


12.9 


166 ± 16 


May 14.03 


VLA 


21.8 


28.5 


133 ± 23 


May 23.96 


VLA 


4.8 


8.7 


191 ± 25 


May 23.96 


VLA 


7.4 


8.7 


248 ± 18 


June 1.06 


VLA 


22.5 


26.5 


118 ±20 


June 12.88 


VLA 


4.8 


13.8 


140 ± 26 


June 12.88 


VLA 


7.8 


13.8 


117 ±17 


June 26.86 


VLA 


4.8 


14.0 


158 ± 19 


June 26.86 


VLA 


7.4 


14.0 


160 ± 15 


July 9.94 


VLA 


4.8 


14.0 


62 ±23 


July 9.94 


VLA 


7.4 


14.0 


98 ±21 


Aug. 21.71 


VLA 


8.5 


37.2 


< 44 


Dec. 16.45 


VLA 


8.4 


25.0 


63 ±20 
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Table 3. PTF GRB Simulation Results 



N 


A^Dct 


A^Ob. 


1 


1311 


11376 


2 


1583 


40101 


3 


228 


5889 


4 


118 


825 


5 


30 


693 


6 


31 


305 


7 


6 


189 


8 


4 


113 


9 


3 


54 


> 10 


26 


426 



